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Slave particle approaches are widely used in studies of exotic quantum phases. A complete 
description beyond mean field also contains dynamical gauge fields, while a simplified procedure 
considers Gutzwiller-projected trial states. We apply this in the context of bosonic models with 
ring exchanges realizing so-called Exciton Bose Liquid (EBL) phase and compare a Gutzwiller wave 
function against an accurate EBL wave function. We solve the parton-gauge theory and show that 
dynamical fluctuations of the spatial gauge fields are necessary for obtaining qualitatively accurate 
EBL description. On the contrary, just the Gutzwiller projection leads to a state with subtle 
differences in the long-wavelength properties, thus suggesting that Gutzwiller wave functions may 
generally fail to capture long-wavelength physics. 
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I. INTRODUCTION 

Over the past decade, the study of spin liquids and 
non-Fermi liquids have been an active theme in con- 
densed matter physics.— A common approach used in 
many of these studies involves the notion of fraction- 
alization where the original particles of a microscopic 
model are substituted with slaves particles coupled to 
a gauge field.— It is often thought that Gutzwiller wave 
functions, constructed by performing projection into the 
physical Hilbert space, are able to capture the correct 
physics. However, it is suspected that such wave func- 
tions may be not sufficient to capture the long wavelength 
properties in important cases with gapless gauge fields, 
e.g. for U(l) spin liquids,—"^"— as the Gutzwiller con- 
struction does not include spatial gauge fluctuations.— 
In a recent study of a hard-core boson model with pure 
ring exchange interactions which we proposed as a can- 
didate model for realizing an Exciton Bose Liquid (EBL) 
phase,—"— we noticed that the EBL can be viewed as 
a special solvable example of a gapless parton-gauge 



system 
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In this work, we shall take up this critical 



issue that Gutzwiller wave functions might not capture 
the spatial gauge fluctuations by explicit demonstrations 
in the EBL context. 

To set the stage for our discussion, we begin with a 
schematic hard-core boson model with ring exchange in- 
teractions which serves the dual-purpose of introducing 
the EBL theory as well as motivating the wave functions 
used in this work. The Hamiltonian defined on the square 
lattice is 



P,rm(r) - bib 



r+mx ^r+mx+ny br+ny, 



(1) 



(2) 



where Pmn(r) are extended ring exchanges on m x n 
plackets and Kmn are amplitudes for these exchanges. 
The Hamiltonian conserves boson number on each row 
and column, and throughout our Hilbert space is the 



sector with equal number of bosons on each row and col- 
umn. We will assume Kmn > 0. Our recent study of 
a model with Kn and i^i2 = K21 found regimes of the 
EBL phasej^i^i while here we are not concerned with 
a detailed realization but rather qualitative aspects, as- 
suming the model Eq. (jlj is deep in the EBL phase. 

In the following, we consider a slave-particle approach 
applied to this problem. Writing each boson operator as 
a product of two parton operators 



bl 



(3) 



we can recover the physical Hilbert space by imposing 
the constraint n(r) = ni(r) = n2(r). We then consider 
states where the hi partons hop only in the x direction 
while 62 partons hop only in the y direction^^ (so a sin- 
gle microscopic boson b indeed cannot hop by itself) , and 
further justify this by noting that the mean-field expec- 
tation value of each ring term in Eq. ([T|) acquires a large, 
negative energy 

(-i^,„„P„„(r))„^f = -/w|Gi(m.T)|'|G2(ny)|'(4) 

Gf,{mfl) = (6t^,6r-HmA,M)m.f. • 

Beyond the mean field slave particle treatment, we in- 
troduce fluctuations into the theory by coupling the two 
parton species to a gauge field a residing on the links of 
the lattice, with opposite gauge charges for the respec- 
tive species. The parton-gauge system is qualitatively 
captured by the following U(l) lattice gauge theory22 

^u(i) - -t5][e^'^-'^-&]:^6r+A,M + H.c.] (6) 
+ hY,4t.~ kY,cos{V ^a\, (7) 

(V-e)r = ^<zX^rM , (8) 

where = ±1 (^ = 1, 2 or x, y) are the gauge charges for 
the partons moving respectively along x and y directions. 
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Equations (jG]) and (O are the respective Hamiltonians 
for the partons and the gauge fields, while Gauss' law in 
Eq. ([8]) imposes a constraint on the physical states. The 
lattice curl and divergence used in Eqs. (O and © are 
defined by 



A,B of Ref. I24 b orrowing from stability analyses of the 
EBL in Refs. Il7l - [l9l) . To study the qualitative effects 
of spatial gauge fluctuations, we consider the following 
parton-gauge Lagrangian which provides a transparent 
starting point for our analysis 



(V X a)r 
(V • e). 



^r-\-x,y Ojry ^r-\-y.x ~r ^rx ; 
^rx ~^ ^rv ~ ^r—x.x ~ ^r— 



flrx , (9) 
r-J/,l/ • (10) 



In the above U(l) gauge theory, the integer- valued 
"electric" field Cr^^ is canonically conjugate to the com- 
pact gauge field on the same lattice link. Dynamical 
fluctuation of these fields arises from the competing terms 
in the gauge field Hamiltonian. In the limit h ^ K,t, 
the electric field vanishes and the Gauss' law reduces to 
ni(r) = n2(r), which projects back into the physical bo- 
son Hilbert space. In this limit, it is possible to eliminate 
the gauge field perturbatively and obtain a Hamiltonian 
for hard-core bosons on the square lattice with ring ex- 
change terms of the type in Eq ([T]), thus establishing 
formal connection between -ffring and -ffu(i)'^ 

As we will argue below, the EBL phase in iJring corre- 
sponds to a "deconfined" phase of i?u(i)i where we can 
ignore the compactness of the gauge field, and treat the 
spatial gauge fiuctuations fully. This is possible in the 
present case due to the powerful bosonization technique 
made applicable by the one-dimensional character of the 
partons and some "dimensional reduction" occuring in 
the systemjiir— i^^TiS^ On the other hand, a different route 
beyond mean field often used in the literature is to ap- 
ply Gutzwiller projection, mostly popular because of its 
numerical tractabilityS£i^ (while gauge theories are of- 
ten intractable). As one can anticipate, this state does 
not know about the spatial gauge field fluctuations and 
fails to reproduce the long-distance properties of the EBL 
phase. 

The paper is organized as follows. In Sec. [Ill we start 
from a Lagrangian formulation and show how the gauge 
theory leads to the EBL field theory, while neglecting 
the spatial gauge field fluctuations leads to a decidedly 
different low energy effective theory. In Sec. Illli we con- 
struct the wave functions used in this paper, and derive 
results for density structure factor and box correlator in 
the harmonic approximation for the wave functions. In 
Sec. lIVl we present our accurate Variational Monte Carlo 
(VMC) calculations for hard-core bosons and show that 
the Gutzwiller wave function indeed realizes a quantum 
state that is distinct from the EBL. In the conclusion, we 
discuss our study more broadly. 



II. EFFECTIVE ACTIONS FOR THE EBL AND 
GUTZWILLER THEORIES 

For the remainder of the paper, we assume a stable 
"deconfined" phase of ifu(i) where we can ignore com- 
pactness of the parton phase variables and compactness 
of the gauge field (stability is discussed in Appendices 



^ = TT [.9"'(5.T^i)' + .9(5x01 - a.)'] + -{d.Oi){dr(j)i) 

ZTT TT 

+ ^ [g~\dye2f+9idyCl>2+ayf] + " (9, ) (^T ^2 ) 

ZTT TT 



(11) 



where the coarse-grained fields (f>^ and 0^ provide a hy- 
drodynamic fiuid description of partons moving in the re- 
spective direction /i, and minimally coupled to the gauge 
field a. The velocity v and dimensionless parameter 
g are convenient parametrization from the bosonization 
literaturej^"— In this formulation, 0^ gives the phase of 
a parton while the dual variable 6^ is related to the par- 
ton density fluctuation through i5n^ = tt^^Q^^^. Wc also 
assume a sizable "stiffness" k for the gauge field (e.g., set 
by the energetics of the boson ring exchanges). Instead 
of introducing the temporal gauge field, we impose the 
following constraint at each lattice site 



(12) 



which allows to recover the physical Hilbert space by 
binding two partons to give the original boson. The con- 
straint is then solved by introducing a field § satisfying 



(13) 



where, as the analysis below shows, -d can be identified 
as the coarse-grained field dual to the boson phase (j) in 
the "bosonization" of the two-dimensional ring exchange 
model in Ref [TtI. 

We first integrate out the fields </)^ and obtain 



-Ccff = —{dxdyd)^ + -^ 
irg Zirvg 



+ -{dr'd){dxay - dyttx) + TridxCiy - dya^f . (14) 

TT Z 

After further integrating out the gauge field a and then 
dropping a less relevant term (dj-'^d)'^, we arrive at the 
following realization of the EBL theory 



C 



1 



EBL 



27r2K 



Tig 



[dxdy-df 



(15) 



where the more general EBL theory is defined by the 
actioE^ii 

§ebl[^?] - i^A^EBL(k,w) |l?(k,L.)|2, (16) 



Xebl(0,w) 7WEBL(k,0) 



\Kky\\ (17) 



for small k^, ky. [Strictly speaking, going from Eq. (|14p 
to Eq. ((TS|) . we need to keep A^EBL(k, w) accurately on 
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the full lines k = (0,fcj,) and {kx,0), i.e., we should not 
drop the naively less relevant term (9t-Vi?)^. However, 
here we focus on long- wavelength effects originating near 
k = (0, 0) and work in a schematic continuum nota- 
tion, while an accurate lattice variant can be found in 
Appendix B of Ref. Hlj . The energy dispersion can be 
obtained from Eq. PT|) and has the form £'k ^ \kxky\. 
This is responsible for interesting properties of the EBL 
phasoii such as specific heat C ^ Tlog(l/r), which 
makes it qualitatively different from sliding or cross- 
sliding Luttinger liquid phases i^"— [Generally, the van- 
ishing of £^k along the lines {0,ky) and {kx,0) can be 
shown to be a consequence of the conservation of bo- 
son number in each row and column of the lattice ring 
model, and is satisfied in this parton-gauge approach by 
construction.] 

Let us now sec what happens if we do not have dynam- 
ical gauge fields. To obtain the resulting Lagrangian, we 
drop the gauge field from Eq. p^ : 



V 

ng 



-Ccutzw = -^[dxdy'dY 



1 



2iTvg 



(18) 

We will view this as a schematic model of what happens 
under Gutzwiller projection, hence the label "Gutzw". 
The corresponding action is 

■^Gutzw 

= ]^^MGnt7.A^,^) |i9(k,w)p,(19) 
XGutzw(k,t^) = —\kxky\^ + —u''\^. (20) 



Try 



nvg 



Here, the energy dispersion is ^ {k^kyl/lkl and the 
distinct behavior in the vicinity of k = leads to low 
energy properties different from the corresponding EBL 
properties. For example, the specific heat vanishes lin- 
early with temperature for the Gutzwiller action, i.e., 
does not have the logarithmic factor log(l/T) found for 
the EBL case. 

The long wavelength properties of the EBL and 
Gutzwiller actions are also different. To give examples of 
other observable consequences, we calculate the density 
structure factor D(k) and box correlator B{x, y) defined 
below: 



B{x,y) 



(21) 
(22) 

(23) 



with (j>r denoting the boson phase variable, 6j ^ e 
Here, (l^kP) can be evaluated for the Gaussian action 
using drir = ir'^dxdy'&r, and we obtain 



for small fc^. ky. The singularity in the structure factor is 
distinct at k = for the two actions. Specifically, at fixed 
ky, D{kx 0,ky) = C{ky)\kx\, with C^BLiky) ^ \ky\ for 
small ky, but CGutzw(fcK) ~ const for small ky. 

Since rir and 4>r are canonically conjugate to each 
other, they satisfy the following ground state minimum 
uncertainty relation 



v/(M^v/O0kF> = l/2, 



(26) 



which allows to obtain the box correlator Eq. ([22]) . Wc 
will focus on the regime |a;| ^ \y\, where we find power 
law decay ^ |a;|^'''-^^ with y-dependent exponents. To 
determine the exponents for all y, we in fact need to have 
details on the (0, ky) line all the way up to the Brillouin 
zone boundary [only the large y limit is determined by 
focusing on the vicinity of k = (0,0)]. For illustrations 
below, we simply take model I?EBL(k) and Z)Gutzw(k) by 
replacing 2|sin(fcy/2)| in Eqs. ([Ml) and For 

the EBL case we fincU^ 



BEBhix,y) 
'7ebl(2/) 



-»;EBL(y) 



1 /27rv r sixC{kyy/2) 



gn Jq sm{ky/2) 



1 2ttv 
— J log(y , 



(27) 
dky (28) 



(29) 



where the last line gives growth behavior for |y| ^ 1. For 
the Gutzwiller box correlator wc find 



BGutzw{x,y) 
jyoutzw 



/2/g 



(30) 
(31) 



which is independent of y in the present Gutzwiller model 
and generally remains finite for any y. 

For finite y and large x, B{x,y) can be viewed as the 
propagator for an exciton of transverse size y. The quali- 
tative difference in the box correlator for large transverse 
size shows that the two actions indeed lead to different 
long wavelength properties. Thus, whether or not one al- 
lows gauge fiuctuation, does lead to effective low energy 
theories with distinct ground state properties. 

We emphasize here that the stability of the EBL phase 
in the ring model given in Eq. ([l} is not the focus of 
this study. Instead, we take the parton-gauge action in 
Eq. (|lip as our starting point and address the question 
whether excluding gauge fiuctuation may lead to a quali- 
tative difference. Note that the Gutzwiller action is only 
a caricature of what happens under the Gutzwiller pro- 
jection and one should use some effective parameter ^off 
rather than bare g. In the next section, we will give 
a more accurate treatment by explicitly constructing a 
Gutzwiller wave function and comparing its properties 
with those of a model EBL wave function. 



DEBL(k) 
£'Gutzw(k) 



2lTV 



kx kf, 



(24) 
(25) 



III. TRIAL WAVE FUNCTIONS 

In this section, we examine the formal properties of 
the Gutzwiller and EBL wave functions and highlight 
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qualitative difFercnces between them. 



A. General Jastrow wave function and harmonic 
approximation 

We first derive expressions for the density structure 
factor and box correlator for a general Jastrow-type wave 
function with a two-body pscudo-potentia l^^'^^ 



^ {{vi}) cx exp 



(32) 



where the indices i,j run over the bosons. In the second- 
quantized notation on the lattice, the wave function can 
be equivalently expressed as 



l^*) cx ^ exp 



\{nr}) 



(33) 

We will shortly see that both the EBL and Gutzwiller 
wave functions have such forms, and their pseudo- 
potentials u{r) will be given later. If we disregard the 
discreteness of the boson number here, we obtain the 
following approximate density structure factor— for an 
arbitrary Gaussian wave function (viewed in variable) 



1 

2uk' 



^u(r)e 



(34) 
(35) 



The box correlator defined in Eq. ((22)) can be calculated 
using 



1 



(36) 



which follows from the boson phase operator (/>r being 
canonically conjugate to the boson number operator n^. 
Again, we have made use of the harmonic approximation, 
that is, we neglect the discreteness of n^, or equivalently 
the compactness of 0r • 



B. EBL wave function 

For a model EBL wave function, we use the pseudo- 
potential from Refs. [2lll22l . which can be motivated by a 
direct "spin-wave" treatment of the ring Hamiltonian in 
Eq. (P, 



1 



Webl e 



zkr 



"EBL(r) ^2 4| sin(fc,/2) sin(fcy/2)r 



(37) 



Note that we exclude lines {kx,0) and {0,ky) from the 
sum. One can also turn this into a convergent integral 



by replacing e^'^^^ by e^''^^ — 1 and e^'^"^ by e*'''"^ — 1; this 
does not change 4'ebl because of fixed particle number 
in each row and column (in our working Hilbert space 
appropriate for the ring models). In principle, VFebl can 
be a smooth function of k but for simplicity here, we take 
it to be a constant. We now use harmonic approximation 
and obtain the following density structure factor and box 
correlator 



^EBL(k) 

Seel (a;,?/) 



Webl 



sm 



sm 



(38) 
(39) 



where Eq. pop holds for large x and fixed y and t^ebl (y) 
is given by 



mBh{y) 



EBL 



2Webl 



" sin'(fc^j//2) 
sm{ky/2) 



dk„ 



2y-l 



HA 



EBL 



log(j/), y > 1. 



(40) 

(41) 
(42) 



Properties Eq. dSH]) with L'EBL(k) \k^ky \ and Eq. ^ 
with rjEBhiy) growing logarithmically with y, are long 
wavelength properties of the EBLJi 



C. Gutzwiller wave function 

To obtain the Gutzwiller wave function, we use the fol- 
lowing wave function for partons confined within a chain 



*cliain {{Xi}) (X exp 



L^2|sin(fc,/2)| 



(43) 
(44) 



(We can again regularize the sum by replacing 6**^=^ by 
gik^x _ 2 since adding a constant to the pseudo-potential 
does not change the wave function for fixed particle num- 
ber in the chain.) This trial wave function has been 
known to capture the energetics as well as Luttinger liq- 
uid exponents of one-dimensional systemsi^^^^ We con- 
struct the Gutzwiller wave function as 



«'Gutzw({rJ) = «'i({rJ) «'2({rJ) 



(45) 



where is the wave function for the 6^ partons confined 
to move within chains oriented in the fi direction. Note 
that Gutzwiller projection has been explicitly imposed in 
Eq. (j45p where both parton species are present at each 
boson location for any given set of {r^}. The Gutzwiller 
wave function indeed has a Jastrow form with the follow- 
ing pseudo-potential 

UGutzw(r-r') = Sy,y uidix-x') + S.j:^x, wid(y-y')- (46) 
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Again disregarding the discreteness of the boson num- 
bers, we obtain the following density structure factor 



D, 



Gutzw 



(k) 



1 



|sin(fc^/2)sin(fc^/2)| 
Wid |sin(fc,/2)| + |sin(fcy/2)| 



(47) 



This has a different singularity at k = compared to 
the density structure factor for the EBL wave function 
in Eq. 

The difference in the structure factors near k = man- 
ifests itself in the fluctuation properties. For a rectangu- 
lar region [0,x) x [0, y), we define the following number 
fluctuation for the total number of bosons in the region 



x—l y—1 



(48) 



The variance of the number fluctuation is readily calcu- 
lated 



L2 



iiin{kxxl2) sin(fey?;/2) 
sm{kx/2) sin(fcj,/2) 



In the limit a; ^ y 3> 1, this has different asymptotic 
forms for the EBL and the Gutzwiller wave functions; 



{SNGutzwix,y)^) 



tt^Webl 



log{x) log(2/), (49) 
(50) 



Equation ()49p shows that such number fluctuation in the 
EBL wave function is strongly suppressed. On the other 
hand, it scales linearly in the region width y for the 
Gutzwiller case while increasing logarithmically with x. 
This reminds of the additivity of variances of statistically 
independent random variables, and the Gutzwiller result 
appears to suggest that, in the absence of gauge fluctu- 
ations, the bosons in adjacent chains are weakly coupled 
compared to those in the EBL phase. 

We now turn to the box correlator and obtain 



Gutzw 



(51) 



for large x and finite y, which is again qualitatively dif- 
ferent from the EBL box correlator in that the exponent 
here does not grow with y. Notice that the Gutzwiller 
result has in fact identical power law to the mean field 
box correlator [see Eq. ([5])] 



ix,y) = |Gi(x£)p \G2{yy)\' 



(52) 



i?/r^"/". (53) 



This therefore suggests that the Gutzwiller projection 
has not provided any improvement over the mean field 
slave particle treatment as far as long-distance proper- 
ties are concerned. [Note that the schematic treatment 
in Sec. HI] leading to results Eqs. ([25]) and ([3T|) might sug- 
gest otherwise if we naively use g = g-ax.i. = """/W^irf there; 



however, such treatment appears to over-emphasize the 
role of the constraint on the long-distance properties of 
the wave functions and we should allow some effective 
Poff instead. We believe the direct approach to the wave 
functions as in this section is more accurate and shows 
that there is no change in the power laws compared to 
the mean field.] 

To conclude our harmonic approximation study of the 
EBL and the Gutzwiller wave functions in this section, 
we have shown that despite the ability of the Gutzwiller 
wave function to realize a quantum liquid, it does not give 
a fully qualitatively accurate representation of the EBL 
phase as defined by Eqs. (fT6|) and (fTT]). In the next sec- 
tion, we enforce hard-core boson condition at each lattice 
site in Variational Monte Carlo calculations and obtain 
numerically exact information for the corresponding wave 
functions, defined in sectors with fixed boson number in 
each row and column. 



IV. EXACT VMC RESULTS 

In this section, we perform exact calculations for the 
hard-core bosons using the wave functions from Sec. IIIII 
We set up Variational Monte Carlo simulations which al- 
low hard-core boson constraint to be imposed exactly. 
We also require fixed boson number in each row and 
column. Since we are only interested in wave functions 
which realize liquid phases, it is important to ensure that 
the variational parameter chosen for each trial wave func- 
tion does not lead to an ordered phase. For concrete- 
ness, we choose density with p = 1/2 (i.e., L/2 bosons in 
each row and each column of L x L lattices), and select 
Webl = 1.5 for the EBL wave function and Wu = 1.5 
for the Gutzwiller wave function so that both wave func- 
tions are deep inside the liquid regimes 4i For the Monte 
Carlo random walks, we allow all possible m x n ring 
moves where bosons hop from occupied sites at r and 
r+mx+ny onto vacant sites at r+mx and r+ny. These 
are the simplest moves that preserve the boson number 
in each row and column and also guarantee ergodicity in 
the Hilbert space of the problem. 

We now present the results of our numerical study. In 
Fig. (U we analyze the density structure factor for each 
of the two wave functions by plotting 



C{ky) 



D{kx , ky) 



sm{kx/2) 



(54) 



kx—l'K / L 



taken at the smallest k^ = 27r/L. This gives a finite- 
size measure of the slope of the density structure factor 
characterizing the V-shaped singularity in the small kx 
limit at fixed ky, and we are further interested in the 
behavior of C{ky) for small ky. In the top panel, we 
obtain the limiting behavior for the EBL wave function 
CEBL(fcy — )• 0) w [fcyl/WEBL, which agrees very well with 
the results derived using the harmonic approximation in 
Eq. ([55)1 (illustrated as a broken line in the figure). 
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FIG. 1: Analysis of the VMC density structure factors at 
long wavelengths for the hard-core boson wave functions. We 
plot C{ky) = D{kx,ky) /\sm(kx evaluated at the smallest 
kx =2ii / L versus ky for system sizes from L = 20 to 80. Top: 
The EBL result shows C{ky) approaching |fcj,|/WEBL for small 
ky. Bottom: The Gutzwiller result shows C{ky) approaching 
the constant l/VKcutzw for any fixed ky ^ upon increasing 
L. The broken lines show the results obtained using harmonic 
approximations in the L ^ oo limit. 



The bottom panel of Fig. [T] shows the corresponding 
analysis of the density structure factor for the Gutzwiller 
wave function. Here, CGutzw(fcy) approaches constant 
l/W'^Gutzw (horizontal broken line) for any finite ky when 
lattice size L ^ oo, which again is in line with the re- 
sult in the harmonic approximation in Eq. (j47p . We also 
examine the ratio of the VMC Gutzwiller density struc- 
ture factor to that in the harmonic approximation (not 
shown), and verify that the VMC data indeed converges 
toward the analytical trend in Eq. (|47)) with increasing 
L. Thus, the density structure factor at long wavelengths 
clearly has a qualitatively different behavior for the EBL 
and Gutzwiller hard-core boson wave functions. 

Figure [2] shows the box correlator B{x,y) versus x for 
y = 1 to 10, measured on a 80 x 80 lattice. The results 
for the EBL and Gutzwiller wave functions are given in 
the top and bottom panels respectively. The data points 



FIG. 2: Box correlator B{x, y) versus x for y = 1 to 10, mea- 
sured on a 80 X 80 lattice. The data are fitted to the ansatz 
in Eq. (|55p . and the results show that the power law expo- 
nent ri{y) increases logarithmically with y for the EBL wave 
function (top), while the exponent is essentially independent 
of y for the Gutzwiller wave function (bottom). 



are plotted together with the best-fit curves using the 
following ansatz 

-viy) 

B{x,y) = Ay 3 sin (^) , (55) 

(56) 

From the two plots shown in Fig. [2l it is clear that 
Eq. ([55)1 provides very good fits for the data. For fixed 
y, the parameter r](y) determines the exponent in the 
power-law relation B{x, y) ^ |a;|~''*^^\ 

For the EBL wave function, the lines fanning out in the 
top panel show that the fitting parameter ri{y) increases 
with y. Anticipating a logarithmic relation from Eq. (j42p . 
we perform an additional data fit to the following 



L 




^■KX\ 




sin 1 










x\- 


-ri{y) 


X 



77(y) = 7log(y), 



(57) 



and obtain 7 — 0.148. This value is very close to 
Webl/tt^ from the harmonic approximation, and there- 
fore suggests that the discreteness of the boson number 
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and the hard-core repulsion in the wave function do not 
significantly alter the long wavelength properties of the 
resulting quantum state when the wave function is well 
inside the liquid regime. 

For the Gutzwiller wave function, the lines running 
parallel to one another in the bottom panel in Fig. El 
show that the exponent of the box correlator is essen- 
tially independent of y. We obtain very good fits using 
ri{y) = 0.477, which is again very close to the correspond- 
ing Gaussian value Wif^/ir. As before, the discreteness of 
the boson number and the hard-core repulsion do not 
modify the long wavelength results of the harmonic ap- 
proximation. 

We have thus shown that the EBL and Gutzwiller 
results obtained using the harmonic approximations in 
Sec, mil remain valid for hard-core bosons here when the 
wave functions are well inside the liquid regimes. But 
more importantly, the exact VMC density structure fac- 
tor and box correlator show that the Gutzwiller projec- 
tion leads to a quantum state that is qualitatively dif- 
ferent from that of the EBL wave function. We discuss 
further implications of this finding in the conclusion. 

V. CONCLUSIONS 

In this work, we compared a Gutzwiller wave func- 
tion with an EBL wave function, both motivated from 
the same parton-gauge action, and found that they re- 
alize quantum liquids with qualitatively different long- 
wavelength properties. This shows that the Gutzwiller 
wave function, which docs not include fluctuations of 
the spatial gauge field, has failed to capture the long- 
wavelength physics. Similar approaches have often 
been used in the studies of quantum spin liquids and 
other strongly correlated systems, and in some cases, 
Gutzwiller wave functions with gapless partons pos- 
sess competitive ground state energies J^ii^i^Sii^"— We 
note that the gauge fluctuations in our case are more 
damped compared to the cases with generic parton Fermi 
surfacesii^"— or Fermi pointsiii^— , and hence are ex- 
pected to be less important than in those cases, but still 
lead to qualitative effects as we have seen. This there- 
fore raises the possibility that Gutzwiller projection may 
generally fail to capture the correct ground state physics. 

Let us also mention some extensions. An interest- 
ing question in the same setting is to compare entan- 
glement in the Gutzwiller and EBL wave functions^i^ 



and examine the effect of including spatial gauge fluctu- 
ations as modeled by the latter. A direct study both in 
the EBL field theory and in the hard core boson model 
realizations^Si^i would be useful. 

In this paper, we used bosonic partons, which we ar- 
gued to be appropriate for the ring models with Kmn > 0; 
in this case the ground state wave function is positive and 
there is no sign problem. On the other hand, for models 
with Kmn < 0, where in general there is a sign problem 
and the ground state wave function has non-trivial signs, 
it appears to be more appropriate to use fermionic par- 
tons, bl = d\^d\^^ We can argue for this either from 
mean field energetics like in Eq. @, or from the con- 
nection between the corresponding II-cvb,^, and i/u(i)'^ 
This construction gives a so-called extremal DLBL state 
from Ref . [13 where fermionic partons form fiat Fermi sur- 
faces in the mean field. A naive bosonization treatment 
of the corresponding parton-gauge system leads to the- 
ory similar to our Eq. pT|) . and hence to an EBL-like 
long-wavelength description. The Gutzwiller wave func- 
tion = '^dx^d-i is also similar to the one in the present 
study, but with specific sign structure from the product of 
the parton Slater determinants. Inspired by the present 
work, we can attempt to crudely account for the gauge 
fiuctuations by replacing the absolute value |^'di^d2l by 
the Jastrow-EBL form while keeping the sign structure. 
Interestingly, while the density correlations are not sen- 
sitive to the sign structure, the boson ring correlations 
are, and in the frustrated case they have faster power 
law decay (for the same density correlations), as can be 
seen already on the mean field level. 

It would be interesting to examine other contexts 
with gapless parton-gauge systems where Gutzwiller- 
type wave functions have been used, such as gap- 
less spin liquids^ '^^i"'^^'^'^!''^^'^^!^^ and more general Bose- 
metals^^i^i^, and see if we can learn how to include 
gauge fluctuations in these cases, even if only on some 
crude level. 
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